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Abstract. Existing (Spitzer Space Telescope) and upcoming (Herschel Space 
Telescope) facilities are deepening our understanding of the role of dust in trac- 
ing the energy budget and chemical evolution of galaxies. The tools we are 
developing while exploring the local Universe will in turn become pivotal in the 
interpretation of the high redshift Universe when near-future facilities (the At- 
acama Large Millimeter Array [ALMA], the Sub-Millimeter Array [SMA], the 
Large Millimeter Telescope [LMT], the James Webb Space Telescope [JWST]), 
and, possibly, farther-future ones, will begin operations. 



Introduction 



Many of the questions left open today by the Cold Dark Matter (CDM) frame- 
work of galaxy formation and evolution are nested into the basic understanding 
of the physical processes underlying star formation, as the baryonic matter sinks 
into the collapsing and merging dark matter haloes. These physical processes 
are the 'rules' for converting gas into stars and for driving the stellar and AGN 
feedback, which affect the evolution of the luminous and non-luminous baryonic 
component of galaxies (this being, ultimately, what we observe). 

The open questions include, among others: the 'missing satellites' problem, 
for which CDM models predict about ten times, at least, more satel lites around 
massive galaxies than what is actually observed ( Moore et al. 19991 ); the 'char- 
acteristic baryonic mass scale' of galaxies, for which the observed luminosity 
function of galaxies follows the Sche chter functional shape and has a charac- 
teristic stellar mass of ~6xl0 10 Mq ( Kauffmann et al. 20031 ). in stark contrast 
with the scale- free, power-law fun ctional shape of the CDM dark halos mass 
function (e.g., Benson et al. 20031 ); the 'angular momentum' (or 'overcooling') 
problem, for which the sizes o f CDM model galaxies are smaller than what ob- 
served (e.g., iGovernatol 120061 ) ; and the existence of bulge-less galaxies, again, 
difficult to produc e within the CDM framework, where bulges are ubiquitous 
(|Maver et al.ll2008l ). 

Within this broad framework, the investigation of the infrared dust emission 
from galaxies offers a window on the questions revolving around star formation 
and feedback that is complementary to the UV-optical, 'direct stellar light' 
approach. Indeed, about 50% of the light observed in the Universe emerges at 
wavelengths longer than a few jum. and is due to dust-reprocessed stellar light 
( Hauser Sz Dwekll200ll ; iDole et all 20061 ). More extreme still is the fraction of 
UV light, a probe of recent star formation, lost to du st absorption, this fraction 
being around 2/3-4/5 of the total UV from galaxies (lCalzettill200ll ). 
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The role of dust emission in tracing recent star formation becomes more 
important as the galaxy luminosity increases, since there is a loose correlation 
between star formation rate (SFR) and the amount of dust extinction mea 



sured in galaxies and galaxy region s (Figure [p e.g ., Wang Sz Heckman] 



1996 



2001 



Heckman et al.1 Il998l: ICalzettil l20Qll: iHopkins et all l200ll : ISullivan et ~ 

Buat et al.l 120021 ; ICalzetti et al.l 120071 ). The most luminous, and most strongly 

star forming, galaxies in the Unive r se do tend to be infrared-bright (e.g . , I Sanders Sz Mirabel 
19961 : ISmail. Ivison. Blainl[l997l : llvison et al]ll998l:lBlain et al.ll2002|1. even af - 



ter taking into account AGN contributions ( Borvs et al.ll2995l : iPope et al.l 2008 
and r evisions to their infrared spectral energy distributions (SEDs; IPope et~ 
2006). At the most basic level, accounting for the fraction of radiation from 



recent star formation absorbed by dust and re-processed in the infrared en- 
ables an accurate census of the total number of stars formed in galaxies across 
cosmic times, and a bett er understanding of the mass assembly of galaxies 
( Hopkins <fe Beacom 2006, and references therein). 

Furthermore, the investigation of the dust components present in the out- 
flows of galaxies, such as the recently detected Poly cyclic A romatic Hydrocar- 



bon ( PAH) emission in the wind of the starburst galaxy M82 (jEngelbracht et al 



2006), may help understand the energetics of stellar feedback, and shed light on 



the 'pollutants' that can survive in the winds and possibly enter the intergalac- 
tic medium. Stellar feedback h as been sugges t ed as the 'culprit' of the mass- 
metallicity relation in galaxies ( Garnett 2002 ; iTremonti et al.l 120041 ) . and the 
dust / gas ratio is correlated with the oxygen abundance in nearby star-forming 
galaxies ( Draine et al.ll2007l ). 



2. Dust Emission and Dust Geometry 



One common application, when measuring star formation in external galaxies, 
is to use the far infrared (FIR) dust emission as a measure of the dust-absorbed 
UV light, which then enables recovering the 'total' UV and quantifying the SFR 
as traced by the emission of young, massive stars. 

Starburst galaxies sho w a well-defined correlation in the FIR/UV-versus- 
UV color plane (Figure l^ lMeurer. Heckman & Calzettil Il999h . The FIR/UV 
ratio is a measure of the UV attenuation suffered by the system: the higher the 
dust attenuation, the larger the amount of UV stellar energy reprocessed into 
the FIR by dust. Variations in the UV color probe the amount of dust reddening 
present in the system. Thus, in s tarbursts, the UV reddening is a good tracer 
of the total UV dust attenuation ( Calzettil l200ll ). 

Quiescent star-forming galaxies and regions show a ~10 times larger spread 
in the FIR/UV ratio than starburst galaxies, at fixed UV color; the spread is 
in the sense that the starburst galaxies for m the upper enve l ope to the quies- 



cently star-forming systems (Figure [2] and Buat et al. _ 20021. 120051: Bell 



Gordon et aD 12004 : iKong et al.l 120041 : ISeibert et al.1 120051 : ICalzetti et al.M20Q5l ). 



2003; 



Thus, a far larger range of UV dust attenuations seems to be present in star- 
forming objects than in starbursts, for the same UV color. This is possible if 
in quiescent star-forming systems the reddening of the UV colors is not only a 
probe of dust, but also of ageing of the individual star-forming regions contribut- 
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Figure 1. The dust attenuation, expressed as the logarithm of the color 
excess E(B— V)=Ay/3.1 (in magnitudes), as a function of the SFR density, 
in units of M Q yr -1 kpc~ 2 , for 164 star-forming regions in 21 nearby galaxies. 
All the galaxies have metallicity close to the solar value. The continuous line 
is the expected trend after combi ning the Schmidt-K ennicutt law with the 
dust-to-gas ratio of the galaxies ([Calzetti et al.| [2007). The dotted lines are 
the 90% boundary to the datapoints. 



ing to the UV emis sion in the system (regions up to an age of ~100-300 Myr, 
Calzetti et all 120051 ). 



An heuristic scenario can be built by recalling that, within starbursts, the 
large number of supernovae can cause the remnants to merge together and be- 
have as a 'collective entity', which can clear, thanks to the large input of mechan- 
ical energy, the interstellar medium within the starburst site and relocate the 
dust to the surrounding areas. In normal star-forming galaxies, such behavior is 
more fragmented and less energetic, as it takes place in individual HII regions, 
and the energy output is not as spatially concentrated as in starbursts. This 
can cause many of the young HII regions to remain dust-enshrouded for longer 
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Log[L x (FUV)/L,(NUV)] 

Figure 2. The FIR/UV ratio (the ratio of the far-infrared to the far-UV 
luminosity) versus the UV color (here expressed as the GALEX far-UV/near- 
UV color) for star burst galaxies and quiescently star-forming regions. The 
starburst galaxies are shown with star symbols. The grey filled triang les are 
star-forming regions within the galaxy NGC5194 (|Calzetti et al.ll2005h . Red- 
der UV colors correspond on average to larger FIR/UV ratios. The continuous 
line shows the best fit to the starburst galaxies, which is well represented by 
a model of a progressively more attenuated (from left to right) constant star- 
forming population. The dotted line shows the same dust attenuation trend 
for a 300 Myr old stellar population; this model represents a reasonable lower 
envelope to the NGC5194 star-forming regions, and to quiescent star-forming 
galaxies in general. 



times than in starbursts. Secular motions will eventually separate them from 
their parent molecular cloud, but the stellar clusters will have aged while migrat- 
ing to areas of lower dust column density. The presence of multiple mechanisms 
for changing the relative geometry of dust and stars in normal star-forming 
galaxies, and the fragmented nature of their star formation (HII regions behave 
like separate entities, and not like a collective one as in starbursts) is likely at 
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the foundation of the observed spread in UV colors for constant FIR/UV ratio. 
Aged clusters can still be bright enough (at least relative to their younger, but 
more dust enshrouded, counterparts) to contribute to the measured UV emis- 
sion in the normal star-forming galaxies, yielding 'red' UV colors (Figure [2]) not 
because of dust attenuation, but because of age. 



3. The Infrared Emission for Tracing Star Formation 

The bolometric infrared emission from galaxies (integrated over the wavelength 
range A ~5-1000 /mi) has been used as a SFR indicator at least since the data 
from the IRAS satellite showed th at dust emission is widespread and signifi- 
cant in galaxies (Soi fer et al. 19861 ). Young star-forming regions are dusty and 



the dust absorption cross-section peaks in the UV, i.e., in the same wavelength 
region where young, massive stars emission also peaks. Although the conneo 
tion between SFR and infrared emission appears a straightforward one (e.g., 
Devereux &: Young 19921 ). cavea ts have been raised for two major simplifications 



that underlie that connection (Hunter et al. 1986; lLonsdale Persson fc Helou 



119871 : iRowan-Robinson Crawford! Il989l : ISauvaee & Thuanl I1992T ). The first 



simplification is about the opacity of a galaxy: not all the luminous energy pro- 
duced by recently formed stars is re-processed by dust in the infrared; in this 
case, the FIR only recovers part of the SFR, and the fraction recovered depends, 
at least partially, on the amount of dust in the system, as well as on geometry. 
Recipes that account for the light of young, massive stars not absorbed by dust 
have been recently proposed ( Kennicutt et al.ll2009l ): these recipes combine the 



infrared luminosity with the observed Ha luminosity, the latter a proxy for un- 
obscured recent star formation, thus producing 'hybrid' SFR indicators. 

The second simplification requires to neglect the heating of the dust by 
evolved, non-star forming populations; however, the latter will also contribute 
to the FIR emission, leading to an overestimate of the true SFR. If more evolved 
populations contribute mainly to the longer wavelength FIR, this extra contri- 
bution may be calibrated, at least for some classes of galaxies. A few approved 
observing projects for the upcoming Herschel Space Telescope plan to investigate 
this problem. 

The connection between infrared emission and SFR also provides the under- 
lying physical mechanism for the correlation bet ween FIR and radio emission in 



galaxies over 5 orders of magnitude in luminosity (jHelou. Soifer. k, Rowan-Robinson 



1985; iYun. Reddv. Sz Condonll200ll ). Recent analyses based on Spitzer data have 



shown that radio images are smoother versions of the infrared (70 /mi) images 
of galaxies, with a correlation length that depends on the the age of the star 
formation event: cosmic rays created in more recent star formati on events have 



not diffuse d significantly in the interstellar medium of galaxies (jMurphv et al 



2006, l2008l ). The fundamental question of how two different processes, one re- 



lated to the heating of dust by stars and the other affecting the propagation of 
the cosmic rays in galaxies, can produce a tight correlation over many orders of 
magnitude remains, however, open. The higher angular resolution images that 
will be provided by the Herschel Space Telescope in the FIR and by the EVLA 
in the radio will hopefully shed light on this fundamental question. 
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Mid Infrared Emission from Galaxies 



Spitzer's high angular resolution (a few arcseconds) in the mid-infrared has 
offered the opportunity to investigate this wavelength range with unprecedented 
detail, thus building on the foundation laid by the Infrared Space Observatory 
(ISO). 

In particular, Spitzer has enabled a detailed investigation of the mid-infrared 
features in emission in the wavelength range ~3— 18 urn fe.g..lSmith et all 2004: 
Engelbracht et al ll2006l:IWu et al.l l2006l: iDale et al.ll2006l : ISmith et al.ll2007l : lDraine et al 
20071 : iGalliano et alJ I2008TJ ). The dust emission in the mid-infrared (A ~5- 



40 /im) range is characterized by both continuum and bands. The continuum is 
due to dust heated by a combination of single-photon and thermal equilibrium 
processes, with the latter becoming more an d more prevalent over the former at 
longer wavelengths (e.g., Draine Sz Li 20071 . and references therein). The mid- 
infrare d bands are generally attributed to Polycyclic Aromatic Hydrocarbons 
(PAH, iLeger &: Pugedfl984 ; ISellgrenlll984l ). large molecules transiently heated 
by single U V and optical photons in the general radiation field of galaxies or 
near B stars (|Li & Draindl2002l : lPeeters. Spoon k Tielensll200l : iMattioda et al. I 
|2005t), and which can be des t royed, fragmente d, or ionized by harsh UV photon 



fields (|Boulanger et al.lll988l : IPetv et al.ll2005l b 



Beyond being interesting in its own merit, the mid-infrared emission from 
galaxies has known renewed interest for its potential use as a SFR indicator 
in deep galaxy surveys, where the Spitzer (and Herschel) far-infrared bands 
correspond to rest-frame mid-infrared emission in high redshift galaxies (e.g., 
Daddi et al1l2005l : iMarcillac et al.ll2006l ; lDaddi et al.ll2007l ). Furthermore, mono- 



chromatic SFR indicators avoid the uncertain extrapolations required by the 
bolometric infrared emission measurement, when only sparsely sampled dust 
SEDs are available (as often is the case for high-redshift data). The physical 
consideration underlying the use of the mid-infrared dust emission as a SFR in- 
dicator is that the dust heated by hot, massive stars can have high temperatures 
and will preferentially emit at short infrared wavelengths. 

For the spectral regions where the dust SED is largely contributed by 
the PAH band emission, the definition of SFR indicators has been less than 
straightforward. The bands have shown to be heated by both recently fo rmed 
dRoussel et al.ll200ll: iFbrster Schreiber et al.l 12004: IDale et aTll2005l. l2007fl and 



evolv ed ( Haas. Klaas Bianchill2002 ; Boselli. Lequeux Gavazzi 2004 ; Bendo et al 
2008) stellar populations. The PAH emission in galaxies shows a strong correla- 
tion with the emiss ion of the cold dust heated by the general (non-s tar-forming) 
stel lar populatio n ( Haas. Klaas Bianchill20"02l : iBendo et al.ll2008l ). In particu- 
lar, iBendo et al.1 ( 2008! ) finds that the 8 fim emission is more closely correlated 
with the 160 [im emission (cold dust) than the 24 /jm emission (warm dust) on 
2 kpc scales in galaxies. Analysis of the Spitzer 8 /im data of very nearby galax- 
ies, like NGC300 and NGC4631, shows that the PAH band emission highlights 
the rims of HII regions and is depressed inside the regions; this suggests that, 
even when closely associated with HII regions, the PAH dust is heated in the 
photo-dissociation regions s urrounding HII regions and like l y des troyed within 
the HII regions themselves ( Helou et al.l 12004 : IBendo et al.l 12006). In general, 
the impact of the non-star-forming populations on the heating of the 8 fim 
emission is at the level of a factor of about 2, in the sense that the 8 /im Spitzer 
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band will be affected roughly b y this amount of uncertainty when used as a SFR 
indicator ( Calzetti et al.l 20071 ) . 

Far more dramatic i s the dependence of the PAH emission on metallicity 
(lEngelbracht et al.ll2005l : IRosenberg et al.ll2006l : IWu et al.l [20061 : iMadden et al 



20061 ) . This dependence implies that regions with values about 1/10 Zq are 



about 10-30 times underluminous at 8 //m relative to regions with solar meta l- 
licity and same SFR (Figured and lCalzetti et al.ll2007l : IRosenberg et al.ll2008h . 
The observed correlation between the PAH 'd eficiency' and the hardness of the 
interstellar radiation field (jMadden et al.ll2006l ) and the prevalence o f 8 um emis- 
sion i n the photo-dissoc iations regions surrounding the HII regions (jHelou et al.1 
20041 ; iBendo et aill2006l ) in nearby galaxies suggests that the strength/hardness 



of the radiation field influences the strength of the PAH emission, possibly by de- 
stroyin g or ionizing the carriers (|Madden et al.ll2006l : IWu et alJ200& iGalliano et aD 
l2008bh . A competing scenario suggests that in low-met allicity galaxies the de- 
ficiency of PAH e mission carriers may be due to de layed formation, rather than 
destruction (e.g.. lDwekll2005l : IGalliano et al.ll2008ah . In this scenario, PAHs are 
formed in the envelopes around carbon-rich Asymptotic Giant Branch stars, 
later injected in the interstellar medium over timescales of a few Gyr, while 
the other dust components are produced from supernova ejecta, which have 
timescales of a few to a few tens of Myr. If low-met allicity galaxies are intrin- 
sically young systems, there will be a delay between the formation of the PAH 
emission carriers and the formation of all other dust components, thus account- 
ing for the observed PAH emission 'deficiency'. However, this scenario counters 
the results from Hubble Space Telescope observations that show that even the 
supposedly least evolved galaxies in our local Universe harbor stellar popula- 
tions that are at least a few Gyr old, a nd have been forming stars for at least 



that long (e.g., Leo A. ICole et al.ll2007l ). Whatever the reason for the observed 



weakness of the PAH band emission in low-metallicity systems, the impact on 
SFR indicators based on spectral regions where these bands dominate the emis- 
sion can be at the order-of-magnitude level, especially if applied to otherwise 
unknown systems (as can be the case in high-redshift surveys). 

Conversely, spectral regions located at longer wavelengths, where thermal 
equilibrium emission starts to dominate over single-photon processes, become 
increasingly more effective at tracing recent star formation. This is the case, 
for instance of the Spitzer 24 fim band, for which many calibrations exist in 
the literature, applicable to HII regions, st arburst-dominated galaxies, and nor- 



mal star-forming galaxies (Figure [H and Calzetti et al. 2005 : Wu et al. 20051: 



Alonso-Herrero et al.ll2006 ± Perez-Gonzalez et al.ll2006l : ICalzetti et al.ll2007l : iRelaho et al 
20071 : IZhu et al.1 l2008l : lRieke et al.ll2009l ). All of these calibrations assume that 



the objects which the SFR indicator is applied to are negligibly contaminated 
by AGN emission. The 24 fim band has a much lower sensitivity to metallicity 
than the 8 ^m emission, decreasing by just a factor 2-4 for a tenfold decrease 
in metallicity. This small decrease is fully accounted for by the increased trans- 
parency (lower dust -to-gas ratio) of the medium for lower metal abundances 
(jCalzettiet "all 120071 ). 

As in the case of the FIR emission, SFRs based on the 24 //m emission 
rely on the assumption that most of the light from massive, young stars is re- 
processed by dust in the infrared. This is generally not the case in most galaxies, 
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Figure 3. The luminosity surface density (LSD=luminosity/area) at 8 /jm 
versus the Pa LSD of 10 starburst gala xies and 220 star-forming regions in 
33 nearby galaxies (|Calzetti et al.| [2007). The 8 /im emission, from Spitzer 
images, has been corrected for stellar contribution. The Pa line emission 
(A=1.876 /im) has been corrected for dust extinction, and is used as an un- 
biased tracer of massive stars SFR. Of the 220 regions, the ~180 regions in 
high-metallicity galaxies, 12+log(0/H)>8.3, are marked in dark grey, and 
the ~40 regions in low-metallicit y galaxies are in light grey. The starburst 
galaxies (Engclbr acht et ail I2005T ) are low-metallicity ones (black squares). 
The continuous line is the best fit to the high mctallicity star-forming re- 
gions (dark grey), with slope 0.94. Models for a young stellar population 
with increasing amount of star formation and dust are shown as a dash line 
(Z=Z(^) and a dot lin e (Z=l/10 Z©), using the stellar populati on models of 
iLeitherer et al.1 (|1999f ) and the dust models of lDraine fc Lil (|2007t ) . The spread 
of the datapoints around the best fit line is well accounted for by a spread in 
the stellar population's age in the range 2-8 Myr. 



and 'hybrid' SFR indicators are likely to be preferable to the use of a mono- 
chromatic SFR indicator. A hybrid SFR indicator that combines the observed 
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Figure 4. The same as Figure [3j with the vertical axis now reporting the 
24 /xm LSD, also from Spitz er. The black asterisks are the Luminous In- 
frared Galaxies (LIRGS) from lAlonso-Herrero etal] (|2006[) . The best fit line 
(continuous line) has slope 1.2. 



Ha recombination line emission and 24 fj,m dust emission can account for both 
the dust-obscured star formation (traced by the dust emission) and the unob- 
scur ed portion of the star formation (traced by the ionizing photons) (Figure \5j 
andlCalzetti et alJl2007l : iKennicutt et al.ll2007l : IZhu et al.l[200§ : IKennicutt et al. 
20091 ). The proportionality factor between the Ha luminosity and the 24 ^m 



emission changes from HII regions/starburst-dominated galaxies to normal star- 
forming galaxies. This is probably a reflection of the fact that the underlying 
stellar population that provides most of the UV-optic al photons (which hea t 
the dust) changes from one type of system to the other ( Kennicutt et al.l 12009). 
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Figure 5. As in Figure El for the combined Ha and 2 4 /nm LSD. The dat a 
follow the 1-to-l line up to the LIRGs, for b/a=0.031 (jCalzetti et alJl2007f ). 
For this SFR indicator, data and models are degenerate in metallicity. Here 
the 24 fim emission traces the dust-obscured SFR, and the Ha emission traces 
the unobscured SFR. 



5. Summary 

ISO and Spitzer have enabled major progress in our understanding of the relation 
between the stellar populations that heat the dust and the various dust emission 
components present in the mid/far-infrared. 

The mid-infrared PAH emission is produced by dust grains transiently 
heated by both recently formed and evolved stellar populations, thus the use 
of this emission to trace SFRs is subject to large (about a factor of 2 or so) 
uncertainties. The same band emission is sensitive to metallicity showing varia- 
tions of factors 10-30 for changes of a factor of 10 in metallicity The application 
of SFR indicators based on those emission bands to distant objects should thus 
be performed with care. 
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As we move to longer wavelengths, thermally heated dust increasingly con- 
tributes to the infrared emission. While still in the mid-infrared region (A <40- 
50 //m), young, massive stars provide most of the heating for the dust, and this 
spectral region can be used reasonably well for tracing SFRs in galaxies. The 
main caveat here is that AGNs also contribute strongly to the heating of the 
dust emitting in the mid-infrared. Furthermore, all infrared-based SFR indica- 
tors need to rely on the assumption that most of the UV-optical light from stars 
is absorbed by dust, an assumption that is reasonably true only for luminous 
(L>10 10 Lq) galaxies. 

The combination of Spitzer 70 fim and radio data has enabled the investi- 
gation of the FIR-radio correlation at higher angular resolution than any earlier 
data, leading to the understanding that radio maps are smoother versions of the 
infrared maps of galaxies, and the 'diffusion' length of the cosmic rays correlates 
with the age of the star forming event. This is an area where the superior an- 
gular resolution of Herschel and EVLA will be able to unravel the underlying 
correlation between the processes that underlie the FIR and radio emission. 

The Herschel Space Telescope, with its superior angular resolution in the 
far-infrared, will also enable the investigation of the role of the different stellar 
populations in the heating of the different dust components that characterize the 
mid/far-infrared emission. These will, in turn, enable accurate determinations 
of SFR indicators using the long-wavelength infrared emission, of dust masses, 
local radiation field intensities, dust opacities and temperature distributions for 
nearby galaxies, and of the grain size distribution and composition. 

The information on the characteristics of the dust emission and underlying 
dust properties that is being and will be, in the near future, obtained for nearby 
galaxies is going to provide the 'toolbox' for the interpretation of distant galaxies 
that will be observed with Herschel, ALMA, SMA, and the USA-Mexico 50- 
meter millimeter telescope, the LMT. 
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